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We study the two-terminal transport properties of a metallic single-walled carbon nanotube with 
good contacts to electrodes, which have recently been shown [W. Liang et al, Nature 441, 665- 
669 (2001)] to conduct ballistically with weak backscattering occurring mainly at the two contacts. 
The measured conductance, as a function of bias and gate voltages, shows an oscillating pattern of 
quantum interference. We show how such patterns can be understood and calculated, taking into 
account Luttinger liquid effects resulting from strong Coulomb interactions in the nanotube. We 
treat back-scattering in the contacts perturbatively and use the Keldysh formalism to treat non- 
equilibrium effects due to the non-zero bias voltage. Going beyond current experiments, we include 
the effects of possible ferromagnetic polarization of the leads to describe spin transport in carbon 
nanotubes. We thereby describe both incoherent spin injection and coherent resonant spin transport 
between the two leads. Spin currents can be produced in both ways, but only the latter allow this 
spin current to be controlled using an external gate. In all cases, the spin currents, charge currents, 
and magnetization of the nanotube exhibit components varying quasiperiodically with bias voltage, 
approximately as a superposition of periodic interference oscillations of spin- and charge-carrying 
"quasiparticles" in the nanotube, each with its own period. The amplitude of the higher-period 
signal is largest in single-mode quantum wires, and is somewhat suppressed in metallic nanotubes 
due to their sub-band degeneracy. 

PACS numbers: 73.63.-b,71.10.Pm,72.25.-b 



I. INTRODUCTION 

Spin transport represents a new branch in mesoscopic 
physics with several technological applicationgiSiMi^, 
e.g. information storage, magnetic sensors and poten- 
tially quantum computation-. While most theoretical 
models are based on Fermi liquid theory, some work has 
been done on strongly correlated ID systems using Lut- 
tinger liquid theorjiSiSiifliiLiS. This work has focused in 
the weak tunneling regime between the ferromagnet and 
the ID system and found that spin transport may pro- 
vide experimental evidence of spin charge separation, one 
of the main predictions of Luttinger liquid theory that 
remains to be observed experimentally in an unambigu- 
ously accepted way. Despite the possible technological 
applications and contributions to the study of spin charge 
separation in strongly correlated systems, very little ex- 
perimental work has been carried out on spin transport 
in ID systemsji^ This work is complicated by the use of 
multi-walled carbon nanotubes, and explored only situ- 
ations with ferromagnetic contacts with parallel or an- 
tiparallel magnetizations. 

Early experimental work with nanotube devices was 
limited by poor contacts between the electrodes and the 
nanotube, and accordingly theoretical models focused in 
the weak tunneling regime. Recently, Liang et al}^ have 
succeeded in fabricating single-walled carbon nanotube 
devices with near-perfect ohmic contacts to the elec- 
trodes. A schematic representation of their experiment 
is presented in Fig.^ These devices are characterized by 
values of the conductance as high as G = 3.7e^//i, close 



to the theoretically predicted higher limilii^ of Ae^/h. At 
temperatures below lOK, the measured conductance ex- 
hibits approximately periodic oscillations as a function 
of the gate voltage. These oscillations are due to Fabry- 
Perot interference - i.e. quantum interference between 
propagating electron waves inside the resonant cavity de- 
fined by the two nanotube-electrode interfaces. In or- 
der to explain their result, Ref. considered a model 
of non-interacting electrons and used the multi-channel 
Landauer-Biittiker formalism to calculate the differen- 
tial conductance as a function of the bias and gate volt- 
ages. They have found qualitative agreement between 
the calculated conductance and their experimental data, 
especially with regard to the variation of the low-bias 
conductance with gate voltage. 

On the other hand, transport experiments on carbon 




FIG. 1: Experimental geometry (from Ref. 13). A single- 
walled carbon nanotube is located on a silicon gate and ox- 
ide layer. The electrodes, which may be ferromagnetic, are 
grown on top of the nanotube. The doped silicon is used as a 
gate electrode to modulate the charge density. The electronic 
transport properties of the nanotube devices were character- 
ized as a function of bias (V) and gate (Vg) voltages. 
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nanotubesiSiiLiSiiS have shown that electrons in nan- 
otubes are strongly correlated and are better described 
by a Luttinger liquid mode P?i^^'^^i^?i^'*i^^ . This implies 
the electrons in these systems do not exhibit Fermi liq- 
uid properties but instead form collective excitations bet- 
ter described by charge and spin like density waves that 
propagate with different velocities. This Luttinger liquid 
behavior drastically changes the charge conductance for 
these systems and it is interesting to know how this af- 
fects the results obtained for the particular setup used 
in Ref. ^3- Furthermore, this setup can be generalized 
to the use of ferromagnetic electrodes, in order to study 
both charge and spin transport in ID electron systems. 

In this paper we investigate the spin and charge trans- 
port properties in ID electron systems with near-perfect 
contacts to ferromagnetic electrodes (the normal metal 
electrodes correspond to the particular case of zero mag- 
netization). We consider both the case of quantum wires, 
i.e. single-channel electron systems, and single-walled 
carbon nanotubes, but mainly focus on the latter one. 
We use a non-interacting Stoner model to treat the fer- 
romagnetic leads and a Luttinger liquid model for the 
nanotube and consider the case of near-perfect contacts 
to the leads, therefore treating backscattering at the con- 
tacts perturbatively. In order to introduce the effect of 
a finite bias voltage, we use the non-equilibrium Kcldysh 
formalism. Following this procedure, we obtain the con- 
ductance, spin and spin current as functions of the gate 
and bias voltages, the external magnetic field and the 
orientation of the magnetization in each lead. We study 
how the strong Coulomb interactions affect these trans- 
port properties and find some features in the Fabry-Perot 
interference patterns that are related to spin charge sep- 
aration. 



II. THE MODEL 

A single-walled carbon nanotube with long-range Cou- 
lomb interactions is well described by a forward-scatter- 
ing model2^^. In this model the Hamiltonian density is 
given by 



a=la=T,i 



Laa 



+ A 



(1) 



where the right/left moving electron operators iliR/Laa 
have the labels a = 1, 2 for the band and a =t, i for the 
spin projection of the electrons in the nanotube and A is 
the interaction strength. The term in the square brack- 
ets corresponds to the electron density. We also consider 
the same problem for a single-channel quantum wire, for 
which there is no sub-band degeneracy and the band in- 
dex can be dropped. Due to the similarities between the 
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FIG. 2: Schematic representation of the model. 
The leads are modeled as ID non-interacting electron sys- 
tems, the Luttinger liquid parameter is therefore g = 1. The 
nanotube (bold line), on the other hand, is described by a ID 
interacting system, in this case g < 1, which corresponds to 
repulsive interactions. The contacts to the leads are modeled 
as two weak backscattering barriers. The two backscatterers 
generate Fabry-Perot interference. 



two cases, we give explicit analytical formulas throughout 
the paper only for the nanotube, but will present results 
for the quantum wire where appropriate. 

The metallic leads are modeled as two semi-infinite Id- 
non-interacting systemsSSiSLS^, which is obtained with a 
position dependent A, constant in the wire and zero in 
the leads. 

We allow for ferromagnetism in the leads using a non- 
interacting Stoner model*'^ (mean-field treatment of the 
magnetization). The Hamiltonian density is TipM = '^o+ 
Km with Ho = HLhi^ = 0) and 



T La 



aaf3 



(2) 

where (Jap are the Pauli matrices and M is the "exchange 
field" , which is proportional to the magnetization. This 
is constant in each ferromagnetic lead, i.e. M \\ rhi for 
X < —L/2 and M \\ rh2 for x > L/2 and in ordinary 
paramagnetic leads AI = 0. In this case, the total system 
corresponding to a nanotube between two ferromagnetic 
leads is described by the Hamiltonian 



H = I dx Hfm 

'\x\>L/2 



dx HhL 

\x\<L/2 



(3) 



The total Hamiltonian H can take a form identical to 
the Hamiltonian in the case of normal metal leads by ap- 
plying the following transformation to the electron field 
operators separately in the left and right leads respec- 
tively 

^R/Lix) ^ C±^/'^--f-/.d-'*^(-')-VK/L(x), X < -I 

^r/l{x) - e±'/'^--fi:/='^-'^^(^')-V^/i(a:), x>^ (4) 

This transformation leaves TYll invariant and TCfm trans- 
forms into Hq. 

We apply the usual bosonization procedure to study 
this modelS^i^. The four electron modes are associated 
to four bosonic modes described by the fields 9aa and 
their duals faa via the bosonization transformation 



^R/La 



1 



V2^ 



(5) 
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where A is a short-distance cutoff. It is convenient to 
consider the following linear combinations of the fields^^ : 
0^,c/s = 75 ± and 6'±,^ = ^ (6li,^ ± 02,^), with 

1 = 1,2 and ^ ~ c,s. This allows us to define the new 
fields 9i = 9+c, which corresponds to the total charge 
mode and is the only interacting mode, and 62 — 
63 = 6-c and 64 = 0-s', with similar transformations for 
the (fi fields. In terms of these new fields the Luttinger 
liquid Hamiltonian density ^ is diagonalized 



where Vp is the Fermi velocity, v is the renormalized ve- 
locity due to the interactions and g is the Luttinger liq- 
uid parameter. In the inhomogeneous model g and v are 
functions of the position: g — 1 and v — Vp in the leads 

and 9 — J v-p+bx/t; ^ ^nd v = Vp / g ui the nanotube. 



27r 



g{dxipif + - {dx9if 
9 



The contacts between the leads and the nanotube are 
(6) modeled by weak backscattering at the contact points, 
the corresponding Hamiltonian density has the form 



2 2 

Wbs = ^ ^ ^ 5{x-Xm) 

m,a,b— 1 a— ±1 



(7) 



where m labels the position of the contacts: Xi/2 — 
(L is the nanotube length) and and v'^l are con- 
stants proportional to the strength of the backscattering, 
= ^'m/ (27rA) and the same for w^. The backscatter- 
ing terms are restricted by symmetry according to charge 
conservation and spin rotational symmetry around the 
axis of magnetization of the ferromagnetic contact. We 
consider only terms of the form V'JjV'l because these are 
the most relevant in the renormalization group sense 
(the scaling dimension in real space of these terms is 
A = ((7-|-3)/4, while the scaling dimension of terms of the 
form is A = 1). Hence if all scattering terms 

are weak, these terms will dominate. It is straightforward 
to extend the present treatment to include the neglected 
interactions, though we do not attempt this here. 

The effect of the magnetization appears only on the 
backscattering term. In the case of near-perfect con- 
tacts to the electrodes, we can treat the backscattering 
Hamiltonian iJbs as a perturbation to the Hamiltonian 
H — i?LL- This procedure is described in section IhH in 
the context of the Keldysh formalism that we use in order 
to account for the effects of the finite bias voltage. 

The gate voltage introduces a term in the Hamilto- 
nian density proportional to pVg — 2/-nVgdx9x, where p 
is the electron density. The constant of proportionality 
relates the voltage applied at the gate with the voltage 
felt by the nanotube and therefore depends on the sam- 
ple. The Hamiltonian density Ti — Ti-hh + 'Hvg becomes 



'M. = TLiA,, after applying the following transformation 
to the 9i field: 9i — > 9i — VgX , where we absorbed a 
constant of proportionality into the definition of Vg for 
simplicity. This transformation needs to be applied to 
the total Hamiltonian, including the backscattering term 
-ffbs, which means that the gate voltage after this trans- 
formation will only contribute to the perturbation Hamil- 
tonian. 

The effect of the external magnetic field is introduced 
via a Zeeman coupling term in the Hamiltonian 



Hh = -h- dx (V'ka^'a/sV'fla/J + V'LaO'a/sV'La/j) ■ (8) 



The contribution of the magnetic field can be transfered 
to the perturbation Hamiltonian iJ^s using a similar pro- 
cedure to the one described above for the gate voltage. 
Taking the z-direction as the direction of the magnetic 
field, i.e. h — hz, the Zeeman Hamiltonian density be- 
comes Hfi = ~{h/TT)dx92 and applying 02 ^ ^2 + Bx 
(with B = 2h/vp), to the Hamiltonian H — Hll + Hh, 
it transforms as H H^l. The results for non-zero 
magnetic field are presented in appendix IbI 
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III. THE NON-EQUILIBRIUM TRANSPORT 
PROBLEM 



Due to the finite bias voltage the distribution in 
this system is not in thermal equilibrium. This non- 
equilibrium situation is studied using the Keldysh for- 
malism (for a review, see Ref. l29l ). To define a non- 
equilibrium initial state, we assume that until some ini- 
tial time to, the system has reached quasi-equilibrium in 



the absence of impurity scattering u\ 



ah 



,,ab 



0. With- 



out impurity scattering, the total number of right- and 
left-moving carriers, Njj, N[^, are separately conserved, 
so that a partial equilibrium can be established with well- 
defined separate chemical potentials for the right and left 
movers. Hence, the system can be described, up to this 
time, by a thermal distribution governed by the grand 
canonical Hamiltonian 



(9) 



where the chemical potentials in each lead are taken to be 
Ml/2 = and A^a/l = /dsn^./^. The right and left 

moving particle densities are given in the bosonization 
procedure by 71^/^ = ^ Y.aa i^'Paa + daa)- Then 



LL 



2V 

dx — dxipi 



(10) 



We emphasize that the appearance of the voltage V in 
Hv does not represent a physical force on the electrons, 
but rather parameterizes their non-equilibrium distribu- 
tion. 

After the initial time to the evolution of the system 
is governed by a different Hamiltonian H, which as de- 
duced in section mis H = iJLL +-ffbs, with TJll given in 
© and Hhs in Q. We expect on physical grounds that 
introducing localized scattering at the ends of the wire or 
nanotube reduces the current, but cannot affect the non- 
equilibrium distribution in the reservoirs. Hence, we be- 
lieve that the prescription of defining the voltage V from 
the initial distribution as is done using (jSJ gives a faithful 
description of ideal leads. According to this prescription, 
a physical observable, represented by an operator O, is 
then calculated from 



(O) 



1 



where 



Tr(e-''^^'e'^*ae-^^' 



-l3Hv 



Trie 



(11) 



(12) 



The difficulty in evaluating such an expectation value is 
that, unlike in a conventional equilibrium calculation, the 
Hamiltonian Hv governing the initial distribution is dif- 
ferent from H, which governs the time evolution. Thus 
such an expectation value cannot be evaluated by equi- 
librium Green's function techniques. 

Instead, we take advantage of the special property of 
Hv that the voltage couples only to Nj^/i^, which are de- 
coupled "zero mode" degrees of freedom. This technique 



has been applied a number of times before to related 
problems22i2ii2^, but to our knowledge the details of its 
derivation have never been published. For completeness, 
pedagogical value, and to highlight the physical assump- 
tions of the method, we include a thorough derivation in 
appendix The correction to (O) due to the backscat- 
tering is given by 



{SO) = ^Tr 



e-I^H^^S\t)dS{t) 



where 



S{t) ^ T exp 



dt'Hj{t') 



(13) 



(14) 



is the evolution operator for a system with the time- 
dependent Hamiltonian Hj{t'). Here Hi{t) is the Hamil- 
tonian in the frame co-moving with the ideal current, 
defined by 



Hi{t) 



Hll + [Hhs]e^. 



i+Vt ' 



(15) 



with V = ^(Nr-Nl). Note that (because [V, iJLL] = 0) 
all the time dependence in Hj (t) is in the backscattering 
term, and is hence easy to handle when working pertur- 
batively in TJbs- 

Eqs. H13I15|I provide a reformulation of the transport 
problem which is particularly convenient for a perturba- 
tive treatment of the backscattering. Note that - because 
the voltage V appears only within _ffbs ~ a direct expan- 
sion of Eq. H13|) in iJbs will involve equilibrium real-time 
propagators calculated with respect to -Sll- We develop 
this perturbation theory using the Keldysh path integral 
formulation. This involves the usual Trotterization of the 
two evolution operators S\ S in Eq. H13|l using coherent- 
state fields denoted 0+ , (p^ for 5* ( "forward branch" ) and 
6^ ,(p~ ( "backward branch" ) for . Further noting that 
Hhh is quadratic and 77bs acts only at the ends of the 
nanotube/ wire, the fields away from x — ±L/2 can be 
integrated out to obtain the Keldysh integral 



(SO) ^ J V[9^{t)ip^{t)]OKe'^'>-'^'^'"^'"'^'^ 



(16) 



with 

i?pcrt = H^s[^+,9+ + 5,iVt]~Hi,s[^-,e-+S,iVt] . (17) 

Here Ok is an appropriate Keldysh representation of the 
operator O, which can be chosen as usual from fields ly- 
ing on either the forward or backward moving branch, 
or any linear combination thereof - see below for con- 
venient choices. The quadratic action 5*0 is a functional 
of 9"^ (t) , (f^^ (t) , which can be determined from the fact 
that it must reproduce the equilibrium correlation and 
response functions for these fields. Indeed, we do not re- 
quire an explicit expression for So, but instead give the 
response and correlation functions, defined by 

C'ix,t;x',t') = {9{x,t)9(x' ,t')) = \{{9{x,t),9{x' ,t')}) , 

R%x,t;x' ,t') = {9{x,t)9{x' ,t')) = 

^-^Q{t-t'){[9{x,t),9{x',t')]) , (18) 
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where we have apphed the standard Keldysh rota- 
tion to the fields 9^ = 6 zL ^6. By construction 

{d{x,t)d{x' ,t')) — 0. The Green's functions involving 
the if fields are defined in a similar way, replacing 9 hy ip 
in the above equations. There are also Green's functions 
that involve both 9 and ip, these are defined by 

C''^ix,t;x',t') = {9{x,tMx',t')) = ^({9ix,t),ipix',t')}) 
R'^'P{x,t;x' X) ^ {9{x,t)<fi{x',t')) = 

^-iQ{t-t')([9{x,t),^{x\t')]) , (19) 

and similar definitions for C"^^ and R^^. Again, by con- 
struction {9{x,t)(f{x',t')) = {(p{x,t)9{x' ,t')) = 0. 

Using the above procedure we obtain (up to additive 
constants that will not contribute to the final result) the 
Green's functions for the 9i fields: 



Rliit) 



(l-a) 



e{t) + y a^fe e(i _ 2kt„) 

a ^ — ' 

k>l 



J(l-a2)^a2fce(i-(2fc+l)<,) , 



2 

1 - a 



k>0 



\ogt' 



E 

k>l 



k>0 



[{2k + l)t,f 



(20) 



where the subscripts label the position of the contacts 
xi^2 (e.g. Cab{t) = C{xa,t;xb,0)) and 



1-.9 
1 + 



and 



L 

V 



(21) 



We also need the Green's functions for the non- 
interacting modes ^2,3,4, and . These are obtained 



from H2U|I by taking a = 


and replacing i„ 


Rnit) = 




<(t) ^ 













(22) 



The Green's functions for the (p fields can be obtained 
from those for the 9 fields given in (|20|l . by replacing g 
by 1/(7, i.e. by replacing a by —a. On the other hand, 
the only ip Green's functions that contribute to the trans- 
port properties studied in the following sessions are those 



that correspond to the non-interacting modes </52,3,4, and 
therefore they are identical to the functions given in H22|l . 

In order to compute the spin transport properties in 
section |V] we also need the following functions for the 02 
and ip2 fields 



R'-'^-{x,t) = sign(x) - Q{t) Q{\x\-v,t) , 
C'-^-(x,t)^^-^ log 



V„t — X 



(23) 



V^t + X 

C"^^^^(x,i) = C"'^'^^(x,t) and i?'^^*'^ (x, t) = R^^'^^{x,t). 

IV. THE DIFFERENTIAL CONDUCTANCE 

In this section we study the charge transport properties 
of ID electron systems and how these are affected by the 
magnetization of the leads and, more importantly, the 
presence of the strong Coulomb interactions. 

We use the procedure described in sectionsHlland llTTI to 
calculate the differential conductance for these systems. 
This is obtained from the expectation value of the current 
in a nanotube, i.e. a four mode ID electron system with 
the Hamiltonian given in as 

I = T.^Aao.^Rao. " V'LaV'L.J = ^i^.^l) ■ (24) 
' TT 

After a lengthy but straightforward calculation we ob- 
tain that the differential conductance G = dl /dV to sec- 
ond order in perturbation theory is given by 



G 



\ rn 



it) 



sm 



Rlm{t) 



with 



cos(W)| 
(25) 



E[«n)' + M'^^ 

mab 

U2 =2 cos{VgL) J2 ["f + I'f Ml ■ M2 



(26) 



and 



(27) 



and similarly for Rab{t). For a quantum wire (i.e. a 
single-channel electron system) these are replaced by 
Cabit) = 2 (Cf((t) + Cfft(i)) and the global normahza- 
tion is divided by a factor of two (since the quantum 
wire has two modes instead of four). 

Eqs. (|25|l and (|26() are valid for zero external magnetic 
field, which is the case considered in this section, the 
equations for non-zero magnetic field are presented in 
appendix 

Equation (|25|l can be easily generalized to arbitrary 
order in perturbation theory, but the time integrals need 
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to be computed numerically. We present the calculated 
conductance to second order for three different physical 
models in Fig. 01 This models correspond to (a) a nano- 
tube with non- interacting electrons, i.e. taking g — 1, 
(b) a quantum wire with g — 0.5 and (c) a nanotube 
with g = 0.25, which is a physically relevant value for 
single-walled carbon nanotubesiS*i2*i2ii2i2^. The effect 
of the interactions is visible in the dependence of the 
conductance with bias voltage, at constant gate voltage. 

The conductance is a quasi-periodic function of the 
bias voltage. At 1^ = 0, for the non-interacting case, 
see Fig. 121(a), this dependence is a cosine function with 
period 2tt /tp. For a quantum wire. Fig. 131(b), there 
are clearly two different "periods" in the oscillations, 
these are related to the two time scales tp = L/vp and 
ty = gtp = L/v. The existence of these two different 
time scales is due to the two bosonic excitations in this 
system: The spin excitation with velocity v-p and the 
charge excitation with velocity w, and is therefore an ef- 
fect of spin charge separation. The same effect appears in 
Fig. 131 (c), but since for the nanotube there are three non- 
interacting modes with velocity Up and only one mode, 
the total charge, with velocity v, it is less visible than 
in the previous example. The most visible effect of the 
interactions in the nanotube, is the enhancement of the 
amplitude in the conductance around V — 0. This ef- 
fect is observable in the experimental data presented in 
Ref.lli 

The calculated conductance at VgL = Tr/2 as a function 
of the bias voltage for a nanotube with different interac- 
tion strengths corresponding to g — 0.25, 0.5 and 1, is 
presented in Fig. 01 It can be seen using Eqs. H25I) and 
(|26|1 that the conductance for this value of the gate volt- 
age only depends on the Green's functions Cn and iin, 
which depends only on ty. As a result, we can clearly see 
in Fig. 21 that the period of the oscillations is 7r/t„, and 
therefore depends strongly on the interaction strength. 
The amplitude of this oscillations is very small except 
for the first oscillation, which is enough to identify this 
effect. 

As for the dependence with the gate voltage, the con- 
ductance is a periodic function, which is modulated by 
Ml • M2, and this is the main effect of the magnetiza- 
tion in the leads on the conductance. In particular if 
J2ui''u2'' < J2^i''^2''^ there is an angle between the two 
magnetizations for which U2 vanishes for any value of the 
gate voltage, in this case the conductance is given also in 
Fig. El 



V. SPIN TRANSPORT 

In this section we study spin transport properties, i.e. 
the spin density in the nanotube and the spin current 
generated by the magnetization in the leads. 

The spin density expectation value in the nanotube, 
calculated using bosonization and the Keldysh perturba- 
tion formalism as described in sections HTIand lTTll is given 



by 



(28) 



aa!3 



For zero magnetic field, it is technically simplest to cal- 
culate from the bosonizcd form 



5" 



1 



(29) 



and then obtain the other two components by rotational 
invariance. For non-zero magnetic field the calculation as 
well as the final results are much more involved, therefore 
and for the sake of clarity we only present the results in 
appendix IbI 
The result is 



+ sm{VgL) ^^U2 Jdt e<^i2(t) gj^ 



iRiiit) 



+U3 /dte'^i"'*' sir 



2^12(0 

cos(yt) 



sm{Vt) 
sm{Vt) 



(30) 



with 



Ul 



E( ah ah iCr ab ab it^r 

\Ui Vi Ml — U2 V2 AI2 



ah 



U2 = Y^ vfvfMi X M2 , 

ab 

J2 (wf M2 + ufvf Ml 



U3 



(31) 



Notice that the spin density does not depend on the 
position in the nanotube, hence the total spin is LS. 

The first term, proportional to ui, is the known non- 
equilibrium spin accumulation effecl!i*2iiii5i&. It is maxi- 
mum for Ml = — M2, when, in the case of identical con- 
tacts, the other terms vanish. This term does not couple 
two backscatterers, is independent of the gate voltage, 
and is an increasing function of the bias voltage. It is 
depicted in Fig.[Sl[a), since it is the only term term that 
corresponds to the component of the spin in the direc- 
tion of Afi a.t Vg ~ 0. The second term corresponds to 
the component of the spin perpendicular to the plane of 
the magnetizations and is depicted in Fig.[3[c), as func- 
tion of the bias and gate voltages. The third term is the 
only one that survives in equilibrium, i.e. at zero bias, it 
is due to the fact that the backscattering strengths de- 
pend on the spins of the incoming and outgoing electron 
relative to the direction of the magnetizations. It is max- 
imum for Ml = M2, when again for identical contacts 
the other terms vanish. This is depicted in Fig. jS^b), at 
9 = 0, 27r. These terms that couple the two backscat- 
terers, and hence depend on the gate voltage, vary with 
bias voltage in a manner approximately described by a 
sum of two periodic functions, with "periods" given by 
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(b). 
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FIG. 3: Calculated conductance for identical contacts and Mi — M2 for (a) a free electron model, i.e. a nanotube with (7 = 1; 
(b) a quantum wire, i.e. a LL with two modes: spin and charge, with g = 0.5; and (c) a nanotube with g — 0.25, as a function 
of bias V and gate Vg voltages (top) and as a function of the bias voltage V at constant gate voltage Vg = (bottom). As can 
be seen, the effect of the interactions is quite appreciable, in particular with the dependence in bias voltage at = 0. The 
"period" of these oscillations is given by 2tv ftp {tp — 2tv in these figures), in agreement with Ref. 113. but in (b) and (c) there 
is another quasi-periodic component in these oscillations, with period given by 2Tr/tv, the presence of these two time scales, tp 
and t„, is a direct result of spin charge separation. 



27: /tp and 2-1: jt^, as discussed in the previous section for 
the conductance. 
The spin current 



is as the spin density calculated from the component 
in its bosonized form (again see result for non-zero mag- 




2 .3 4 

V (bias voltage) 



FIG. 4: At constant gate voltage: VgL — 7r/2, the period of 
oscillations is n/gtL, i.e. depends strongly on the interaction 
strength. In order for this effect to be clearly visible, we scaled 
and shifted the functions differently, therefore the values of 
the G axes are not meaningful. 



netic field in appendix IB|) 



TT 



(33) 



It is not well-defined at the contact points because the 
backscattering term in the Hamiltonian Q does not con- 
serve spin, and therefore it has different expressions in the 
nanotube and the leads. 

The spin current in the left (-I-) and right (— ) leads is 
given by 



Js =Ui ydte<^"(*) sin 



sin(W) 



[± sin[VgL) U2 + cos{VgL) u^] 



sm{Vt) (34) 



and in the nanotube by: 



+ cos{VgL) \ U2 /die'^i^(*'si 



sin(W) 
1 



Rl2{t) 



(35) 



Us 



die'^i^^*^ sin 



Rl2{t) 



sm{Vt) 
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FIG. 5: Calculated spin for a nanotube with g — 0.25, the component of the spin in the direction of the magnetization Mi as 
function of the angle between the two magnetizations 6 and the bias voltage V, at constant gate voltage (a) = and (b) 
Vg — 7r/2; and (c) the component of the spin perpendicular to the plane formed by the two magnetizations as a function of the 
gate and bias voltages. 



with U2 and defined in (|31() and 

^4 = E<>™^^™- (36) 

mab 

Similarly to the results for the spin discussed above, 
the first term, which only involves one backscatterer and 
is independent of the gate voltage, corresponds to the 
usual spin injection effect. It is an increasing function of 
the gate voltage and is maximum for Mi = M2. This 
can be seen in Fig. IH^b), since it is the only term that 
does not vanish in the direction of Mi at VgL — Tr/2. At 
Vg — 0, the terms proportional to ui and contribute 
equally to the component of the current in the direction 
of Ml, the result for this case is presented in Fig. El^a). 
The second term, proportional to Mi x M2 corresponds 
to an exchange interaction between the magnetizations 
of the leads, mediated by the nanotube. It has opposite 
signs in the two leads and it is shown in Fig.|H|[c). 

VI. CONCLUSIONS 

We studied the charge and spin transport properties of 
ID systems, e.g. quantum wires and carbon nanotubes, 
focusing on the latter. We considered the case of nearly 
perfect ohmic contact between the ID system and the 
electrodes and included the strong Coulomb interaction 
via a Luttinger liquid model. We found important ef- 
fects on the transport properties of these systems that 
are due to the Coulomb interactions. These appear in 
the dependence with bias voltage. In particular, the con- 
ductance is enhanced at low bias voltage, furthermore it 
is an oscillatory function where we can distinguish two 
quasi-periodic components, with periods that are related 
to the two velocities of the excitations of a Luttinger liq- 
uid, V and Up. This effect is therefore a direct consequence 
of spin and charge separation. It is clearly visible in sin- 
gle band quantum wires. In nanotubes, the amplitude of 



the higher period component is reduced by the presence 
of three (as opposed to one) neutral modes. Still, we can 
find evidence of the two velocities v and Vp by compar- 
ing the dependence of the conductance with bias voltage 
for two different gate voltages {VgL — and tt/2). It is 
perhaps worth noting that, for the case of non-magnetic 
leads with symmetric contacts, the conductance formula 
involves only two unknown parameters: the overall am- 
plitude of the backscattered current, and the Luttinger 
parameter g, both of which can be simply estimated. 
Nevertheless, a non-trivial functional dependence upon 
bias voltage is predicted. 

The spin and spin current have one component in the 
plane of the magnetization that does not couple the two 
leads and is therefore independent of the gate voltage. 
This term should be understood as arising from incoher- 
ent spin injection at each contact. It is a monotonic func- 
tion of the bias voltage, and corresponds to the known 
spin accumulation (for the spin) and spin injection (for 
the spin current) effects. The other components that 
couple the two leads, and therefore depend on the gate 
voltage, are backscattering processes occurring with co- 
herence between the two contacts. These oscillate with 
the bias voltage, in a manner approximately described as 
a sum of two periodic components, with periods related 
to the two velocities of the excitation of the Luttinger 
liquid. The amplitude of the higher-period component is 
largest in a single-channel quantum wire, and somewhat 
suppressed in nanotubes by the sub-band degeneracy. 
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APPENDIX A: DERIVATION OF EQ. llTSt 

In this appendix, we derive Eq. (|13() . In particular, 
we consider a large periodic system of size L, where ulti- 
mately L — !■ cx). We define the right /left-moving combi- 
nations 



•/"ifl/L = ± I 



(Al) 



In the system of size L we can decompose into finite 
wavevector and "zero mode" components. In particular, 
for the total charge fields, we define 



2t:Nrx 
[ 



(A2) 
(A3) 



where 4>iR/Li^) contains the non-zero momentum modes 
of the (piR/L fields. With these definitions, the zero-mode 
variables form two canonically conjugate pairs: 



[Nr, $l] 



(A4) 



[Nr,Nl] = [^r,^l]=0 



Moreover, Nr/]^,^r/]^ commute with <piR/Lix) and all 
fields associated with channels 2,3,4. 

Since the interactions which transform the system from 
a Fermi liquid into a Luttinger liquid (Eq. 1^) exist only 
for |a;| < L/2, they do not affect the zero-mode terms in 
the Hamiltonian. Hence one may separate 



V 



(Nr^Nl) , (A5) 



where -ffLL is the Luttinger liquid Hamiltonian, Eq. lO, 
with the zero-mode terms subtracted, i.e. with ipi — s- (pi 



and 01 — > 6*1. We then see, using the independence of the 
zero mode variables, that the unitary operator 



Uv 



(A6) 



generates the following transformation Nr/j^ ^r/l i 
V\-/{A'Kv^), hence 



V ' 



(A7) 

where C = [3\-/{^t:v^) is an unimportant constant. In- 
serting this into Eqs. 111I12II . one obtains 

{O) ^ 



with 



and 



Zi,T. = Tr ie-P"- 



V=^{Nr^Nl) 



(A8) 
(A9) 

(AlO) 

(All) 
(A12) 
(A13) 



For the operators of interest, 

UliUv = lo + I , 

UlSUv = S , 

uIJsUv = Js , 

where Iq — A{e'^/h)V — (2/7r)T^ is the current which 
would flow in an ideal nanotube in the absence of 
backscattering. Defining SO = O — Oq, with Oq — Iq 
ior O — I and ©0 = for O = S, Jg, one has then 

(^0) = ^Trfe-'='«--e*(^+^)*ae''(^+^'*) . (A14) 



We then apply the formula 



-it{H+V) 



-itV 



T exp 



dt'Hi{t') 



(A15) 



to arrive at Eqs. (|13|) ff. given in the main text. 
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APPENDIX B: DIFFERENTIAL 
CONDUCTANCE, SPIN AND SPIN CURRENT 
AT NON-ZERO MAGNETIC FIELD 



The differential conductance including the magnetic field is still given by Eqs. (|25|l and H26|) with only the following 
change in U2 



where h = h/h. 



U2 ^2cos{VgL) |uf + Ml • M2 + {Mi ■ h){M2 ■ h) {cos BL - 1) } 

ab 

+ sm{BL) sm{VgL) ^ [wf (AI2 • h) + ufvf{Mi ■ h) 



(Bl) 



The total spin of the nanotube is 



S ^-hBL - — ^(Uufvf'Mi - ufvfM^ i sinBL + (uf'vf{Mi ■ h)h - ufyfi^h ■ h)hj {L ~ ^ sin BL) 

^ ab 

+ (ul\f{Mi xh)- ufvf{M2 X /i)) i (1 - cosBL)| I e*^'^ sin {\Rii) sin{Vt) 



+ 



(wf M2 + ufwf All) ^ sinBL + (uf (M2 • h)h + ufvf'{Mi ■ 'h)h^ {L cos BL - ^ sin BL) 



sinl^L 



ufufhL sin BL cosV„L - vfv"-'' 



+ 



{Ml ■ h)M2 + {M2 ■ h)Mij -g (1 - cosBL) 
+ (Afi • h){M2 ■ h)h {LsinBL- f (1 - cosBL)) cosVgL^ J^e*^'^ sin(iili2) cos(V<) 

(uf (M2 xh)~ ufvf{Mi X /i)) -|(1 - cosBL) cosT^i + • h){M2 x h) - (M2 ■ h){Mi x h) 

X {L- ^ sinBL) + Ml x M2L] sinVgLl [ e'^^^ sin (5-R12) sin(yt) ) . (B2) 
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The spin current in the nanotube (Eq. IB3|) and the leads (Eq. IB4|) is 
Js = Y^luf'vf'Mi cosB{x - xi) + ufvfli'h cosB{x - X2) 

ab ^ 

+ ufvf{Mi ■ h)h [1 - cosB{x ~ xi)] + uf (M2 • h)h [1 - cosB{x - X2)] 

+ ufvf{Mi X h)smB{x~xi)~ufv2''{M2 x /i) sinB(x - X2)| ^e*^" sin (iiln) sin(V^t) 

+ E| [<''«2^(^2 • h)h [cos BL - cos B{x ~ xi)] + ufvf{Mi ■ h)h [cos BL - cos B{x - X2)] 

+ ufvfM2 cos B{x - xi) + Ml cos B{x - X2)\ cos VgL 

(Ml • k){M2 ■ h)h [sinBL - smB{x - xi) + smB{x - X2)] 
+ (Ml • ^)M2sinB(a; - xi) - (M2 ■ ^)Mi sinB{x ~ X2) 



ab 



+ vfvf 



+ uf ufh sin BL 



''1 "2 



sin l^Lj J e^^-'- sin {\Ri2) sin(FO 
E| [^if^'f (A^2 X h) smB{x - xi) - u^''vf{Mi x h) smB{x - X2) 

(Ml • h){M2 xh)[l- cosB{x - xi)] - {M2 ■ h){Mi xh)[l- cosB{x - X2)] 



sin VgL 



ab^ab 



+Mi X M2 



cosl^i I / e*^!^ sin (5-R12) cos(V^t) , 



(B3) 



Js = J2lufvfMi cosB{x - xi) + ufvfM2C0sB{x - X2) 



+ uf (Ml • [1 - cosS(a; - xi)] + ufvf{M2 ■h)h[l- cosB{x - X2)] 

T [ufwf (Ml X /i)sinB(a;-xi)+'ufi;f (AI2 x h)iiinB{x - X2) | ^e'^" sin (ifin) sin(V^t) 

E{ ["f^fC^z • h)h [cos - cos B{x ~xi)\+ uf (Afi • h)h [cos BL - cos B{x - X2)\ 



cosl^L 



ufvfM2 cosB{x - xi) + ufvfMi cosB{x - X2) 

(Ml • h){M2 ■ h)h [sinSL - sinB(a; - xi) + s\nB{x - X2)] 



+ (Ml • h)M2 sinB(x - xi) - (M2 • h)Mi su\B(x ~ X2) 



sin VgL 



cosl^L 



uf (A/2 X h) smB{x - xi) + <''(Afi x h) smB{x - X2) 

{Ml ■ h)(M2 xh)[l- cosB{x - xi)] - {M2 ■ h){Mi xh)[l- cosB{x - X2)] 



T vfvf 



+Mi X M2 



s\nVgL + ufufhsinBLsinVgL \ / e<^i^ sin (ii2i2) sin(yt) , 



(B4) 



where the =F sign correspond to lead 1,2 respectively. 

One can verify that the spin current and spin density are not independent and are in fact related by the precessional 
equation of motion, 



dtS + d^Js = -2h X S. 



(B5) 



In the steady state, {dtS) — 0, so one has 



9:.(J.(X)) =-2/lX (Six)). 



(B6) 
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